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I.  INTRODUCTION 


The  performance  of  an  LMS  (least  mean  square)  adaptive  array  [11 
can  be  influenced  by  the  modulation  on  an  interference  signal.  For 
example,  pulsed  interference  can  make  the  weights  in  the  array  alternate 
back  and  forth  between  two  sets  of  values,  one  when  the  interference  is 
on  and  the  other  when  it  is  off  [2],  Interference  modulation  has  two 
deleterious  effects  on  the  array.  First,  it  can  make  the  array  modulate 
the  desired  signal.  Second,  it  can  cause  the  array  output  SINR  (signal  - 
to-interference-plus-noise  ratio)  to  vary  with  time.  In  a  digital 
communications  system,  such  SINR  variation  usually  results  in  an 
Increased  bit  error  probability.  For  these  effects  to  occur,  the 
modulation  parameters  must  be  properly  chosen.  For  example,  the 
modulation  rate  must  be  commensurate  with  the  array  speed  of  response. 

In  previous  studies,  the  authors  have  examined  the  effects  on  an 
adaptive  array  of  interference  with  pulsed  modulation  [21,  with 
sinusoidal,  suppressed  carrier  envelope  modulation  [31  and  with 
arbitrary  periodic  envelope  modulation  [4],  Each  of  these  studies 
involved  interference  with  envelope  modulation  but  not  phase 
modulation. 

The  purpose  of  the  present  report  is  to  extend  the  earlier  work  to 
handle  interference  with  phase  modulation.  We  assume  an  adaptive  array 
receives  an  interference  signal  with  a  periodic  but  otherwise  arbitrary 
phase  modulation.  We  present  a  method  of  solving  for  the  resulting 
time-varying  array  weights.  Using  this  method,  we  then  study  as  an 
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example  the  behavior  of  a  2-element  array  that  receives  an  interference 
signal  with  sinusoidal  phase  modulation.  We  show  how  each  interference 
signal  parameter  affects  the  desired  signal  modulation,  the  SINR 
variation  and  the  bit  error  probability. 

In  Section  II,  we  formulate  the  problem  and  present  the 
mathematical  technique  for  solving  for  the  array  weights.  In  Section 
III,  we  describe  the  2-element  array  and  the  phase  modulated 
interference  signal  and  present  our  numerical  results.  Section  IV 
contains  the  conclusions. 

II.  FORMULATION  OF  THE  PROBLEM 

Assume  an  adaptive  array  consists  of  J  isotropic  elements  with  half 
wavelength  spacing,  as  shown  in  Figure  1.  Let  xj (t )  be  the  analytic 
signal  received  on  element  j.  xj(t)  is  multiplied  by  a  complex  weight 
wj  and  then  summed  to  produce  the  array  output  s(t).  The  array  weights 
are  controlled  by  LMS  (least  mean  square)  feedback  loops  [1],  which 
obtain  each  weight  wj  by  integrating  the  product  of  xj(t)  with  the  error 
signal  e(t).  The  error  signal  is  the  difference  between  the  reference 
signal  r(t)  and  the  array  output  s(t).  The  array  weights  satisfy  the 
differential  equation 

+  k*W  =  kS  ,  (1) 

Ht 

where  W  =  [w^,  W2,  ...,  wj]T  is  the  weight  vector,  t  is  time,  *  is  the 
the  covariance  matrix. 
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$  =  E(xV)  . 


(2) 


S  is  the  reference  correlation  vector, 
S  =  E[x*r (t) ]  , 


(3) 


and  k  is  the  IMS  loop  gain.  In  these  equations,  X  is  the  signal  vector, 
X  -  [x^(t),  X£(t ) ,  ...,  Xj(t)  ]  t  (4) 

T  denotes  transpose,  *  complex  conjugate  and  E(*)  expectation. 

We  assume  that  a  desired  signal  and  an  interference  signal  are 
incident  on  the  array  and  that  thermal  noise  is  present  in  each  element 
signal.  The  signal  vector  then  contains  three  terms. 


X  =  X.  +  X.  +  X 
din 


(5) 


where  Xq,  Xi ,  and  Xn  are  the  desired,  interference  and  noise  vectors, 
respecti vely . 

Let  the  desired  signal  be  a  CW  (single  frequency)  signal  incident 
from  angle  0q  relative  to  broadside.  (6  is  defined  in  Figure  1.)  The 
desired  signal  vector  is  then 


jKt'M'd).. 

Ad  e  Ud 


(6) 


where  Aq  is  the  amplitude,  wg  is  the  carrier  frequency,  1>q  Is  the 
carrier  phase  angle,  and  Uq  Is  a  vector  containing  the  Interelement 
phase  shifts, 
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>  •  •  •  » 


(7) 


=  Cl,  e 


-j  4>d 


-jK*d  _T 
e  J 


where 


*d  =  *  sin  9d  *  (8) 

and  where,  to  simply  notation,  we  let 

K  =  J  -  1  (9) 

We  assume  ^  is  a  random  variable  uniformly  distributed  on  (0,2m). 

Next,  assume  an  angle  modulated  interference  signal  arrives  from 
angle  0-j .  The  interference  signal  vector  is 


X,  =A,eJK*  +  +,)Ul(t) 


'1  T 


with 


U.(t)  = 


/ 


V 


jTi(t) 


jCYi  (t -2T -j )  -24>-f  ] 


\ 


j[Ti  (t— KT-f  )-K4>0 


J 


(10) 


(ID 
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where  -y-j  ( t )  is  the  angle  modulation  as  received  on  element  1,  A-j  is  the 
interference  amplitude  and  ^i  is  the  carrier  phase  angle.  T^  and  are 
the  interelement  time  delay  and  carrier  phase  shift, 

T.  =  Jl_  sin  9.  »  (12) 

1  “>o  1 


and 


=  u  T^  =  iTsi  n  6_j 


(13) 


We  assume  il'i  is  a  random  variable,  uniformly  distributed  on  (0,2n)  and 
statistically  independent  of 

Finally,  we  assume  each  element  signal  contains  a  zero-mean, 
independent  Gaussian  thermal  voltage  n j (t )  of  power  cj2.  The  noise 
vector  is 

~  T 

Xn  =  ^nl(t)»  "2^*  •••»  nj(t)]  ,  (14) 


where 


E[nj  (t)nk(t)]  =  a26jk  ,  1<  j,  k  <  J  t  (15) 

with  6j|c  the  Kronecker  delta.  The  nj(t)  are  assumed  statistically 
independent  of  both  tly  and  <|»-j . 

Under  these  assumptions,  the  covariance  matrix  in  (2)  is  the  sum  of 
a  desired,  an  interference  and  a  noise  term, 

$  =  +  <t-j  +  1>n  .  (16) 
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The  desired  signal  term  is 
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To  compute  the  reference  correlation  vector  S  in  (3),  we  must 
specify  the  reference  signal  r(t).  In  practical  applications,  the 
reference  signal  is  derived  from  the  array  output  [5,6,71.  To  make  the 
array  perform  properly,  it  must  be  a  signal  correlated  with  the  desired 
signal  and  uncorrelated  with  the  interference.  Here  we  assume  the 
reference  signal  to  be  a  replica  of  the  desired  signal. 


r(t)  =  Ap  e 


j(u0t+td) 


(21) 


where  Ar  is  its  amplitude.  Equation  (3)  then  yields 


S  =  A  A  .  IJ . 
r  d  d 


(22) 


When  (16)  and  (22)  are  substituted  into  (1),  we  find  that  the 
weights  satisfy  a  system  of  differential  equations  with  constants  on  the 
right  but  with  time-varying  coefficients  (due  to  3>i).  In  the  rest  of 
this  section,  we  develop  a  method  for  solving  this  system  of  equations. 

To  solve  (1),  we  do  several  things.  First,  we  write  the 
interference  covariance  matrix  in  a  more  compact  form.  We  define  J 
vectors, 

Vq  =  (1,  0,  0  ...  0,0)T 

Vj  =  (0,  e"j<\  0,  ...  0,0)T  , 

-  j  2  4>i  T 

V2  =  (0,  0,  e  \  ....  0,0)'  , 
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-jK<f>i  T 

V„  =  (0,  0,  0  ...  0,  e  ) 


Note  that  these  J  vectors  form  an  orthonormal  set. 


V  V.  =  6..  ,  0  <  j,k  <  K  . 

j  k  jk 


By  using  these  vectors,  $-,•  in  (18)  may  be  written 


4>.  =  A?  I  l  V*  VT  . 

i  1  p=0  q=0  P  P 


so  (1)  becomes 


K  K 


+  k  [A2  U*  uj  +  A2  I  l  ejxPq(t)  V*  vj  +  a2l]W  =  kA  A  U* 
L  d  d  d  1  p=0  q=0  P  q  r  d  d 


Next,  we  normalize  (26).  Dividing  by  k a?  gives 


K  K 

dW(t')  +  [•  +  e.  I  l  eJXpqft')  V*  VT]Wft'l  =  Ar  /FT  U*  - 

TF  [  2  i  p=0  q=0  P  qJ  a  d  d 


where 


*2  =  1  +  Cd  Ud  Ud  ’ 


r 


^  =  input  si gnal -to-noi se  ratio  (SNR)  per  element, 

0 


£.  =  ^  =  input  interference-to-noise  ratio  (INR)  per  element, 

a 2 


and  where 

t'  =  ka^t  =  normalized  time  .  (29) 

a 

Also,  we  note  that  the  constant  _L  on  the  right  will  just  appear  as  a 

a 

scale  factor  in  the  solution  for  W(t').  It  has  no  effect  on  the  array 

a 

output  SI  NR  to  be  discussed  below.  Hence  we  arbitrarily  set  _L  =  1  to 

a 

eliminate  it.  Equation  (27)  is  then 


dW(t')  + 
“dT 


[*5 


*  h 


l  I  eJ  Xpq^') 

p=0  q=0 


\  vJ]w(t')  =  K  U 


* 

d  wd 


(30) 


Next,  in  order  to  be  able  to  solve  (30),  we  make  the  important 
assumption  that  the  angle  modulation  Yj(t)  is  periodic.  If 

j  Xpq  (t 1 ) 

Y-j  (t )  is  periodic,  the  functions  e  are  also  periodic. 

Therefore  we  may  expand  each  of  them  in  a  Fourier  series.  In  order  to 
simplify  notation  later,  we  shall  also  include  the  constant  £i  in  this 
expansion.  Thus,  we  write 


j  Xpq ( t 1 ) 

f>i  e 


l 

£  =  -“ 


jdo)  t 


'tpq 


(31) 
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and  w  Is  the 
m 


...  jXpq(t') 

where  f^pq  is  the  Fourier  coefficient  of  Sf  e 
normalized  fundamental  frequency  of  Yi(t'),  i.e., 

«)•  =  “V  2  •  (32) 

m  'ka 


where  is  the  fundamental  frequency  of  Yi(t).  We  note  for  later  use 
that  if  p=q,  then  Xpq(t')  =  0,  so  the  series  in  (31)  contains  only  a 
zero  frequency  term,  i.e., 

Vpp  ~  V  ‘Vo  •  (33) 

Also,  it  is  easily  shown  that 


(34) 


and 


V(p+l)(q+l)  "  Vpq 


I  I 


(35) 


In  addition  to  assuming  Y-j(t)  periodic,  we  shall  also  assume  that 
the  bandwidth  of  the  interference  is  finite,  i.e.,  we  assume  that  only  a 
finite  number  of  terms  are  nonzero  in  (31).  Suppose  the  coefficients  in 
(31)  are  zero  for  |fl|  >  L,  where  L  is  some  integer. 

Then 

jXpq(t')  L  jnut't ' 

Si  e  =  l  f*pq  e  *  .  (36) 

*=- L 
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Since  (30)  is  a  linear  differential  equation  with  periodic 
coefficients,  the  solution  for  W(t')  will  be  a  periodic  function  of  time 
after  any  initial  transients  have  died  out  [81.  In  this  paper,  we  shall 
not  attempt  to  solve  for  initial  transients,  but  shall  consider  only  the 
periodic  steady-state  solution.  Once  any  initial  transients  are  over, 
W(t')  can  be  written  as  a  Fourier  series, 

*  inw't ' 

W(t')  =  }.  Cn  e  m  ,  (37) 

n=-» 


where  Cn  is  a  vector  Fourier  coefficient.  Substituting  this  series  into 
(30)  and  enforcing  the  resulting  equation  for  each  frequency  component 
separately  gives 


(Vjnu)mI)Cn 


L  K  K 

+  I  l  l  f  V*  VT  c  = 
*«-L  p=0  q=0  *P9  P  0  n-i 


★ 

?,U  .6  ,  -"  <  n  < 

d  d  no* 


(38) 


Equation  (38)  is  an  infinite  system  of  vector  equations,  one  for  each  n. 
To  solve  for  the  Cn,  we  first  assume  that  there  is  some  integer  N  such 
that  the  Fourier  coefficients  Cn  are  negligible  for  |n|  >  N.  In  other 
words,  we  assume  W(t')  can  be  adequately  approximated  by  a  finite  sum 


W(f) 


N 


l  ^n 
n=-N 


jno)  t 


(39) 


Such  an  approximation  is  reasonable  because  the  feedback  loops 
controlling  the  array  weights  in  W(t')  are  lowpass  filter  loops  that 
cannot  respond  above  a  certain  speed. 
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\ 


I 
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If  we  set  Cn=0  for  |n|  >  N,  then  (38)  yields  a  finite  system  of 
equations  for  the  remaining  Cn.  Each  vecton  Cn  has  J  scalar  components, 
so  the  result  is  a  system  of  (2N+1)J  linear  equations  in  the  unknown 
scalar  components.  One  could  solve  for  W(t‘)  by  solving  this  system  of 
equations  numerically. 

However,  rather  than  work  with  (38)  directly,  we  shall  instead 
transform  variables  first.  We  do  this  because  the  resulting  equations 
have  a  more  systematic  form.  Specifically,  we  shall  express  each  Cn  in 
terms  of  its  components  along  the  vectors  in  (23).  (Since  each  Cn 
has  J  components,  the  J  orthonormal  vectors  can  serve  as  basis 
vectors.)  We  write 


C 


n 


(40) 


where  the  are  scalar  coefficients,  cy,^  is  the  component  of  Cn 
along  the  unit  vector  .  Substituting  (40)  into  (38),  multiplying 
the  result  on  the  left  by  (for  a=0,l,2,  ....  K)  and  using  (24)  yields 

d 

the  equation 


L  K 


inw  a  +  I  Q  .  a  .  +  I  £  f.  a  „  =  /iT  (V^  U  .)6  , 

J  m  n,a  |<=o  ak  n,k  **_{_  q»o  *aq  n-i,q  d  11  a  d'  no  * 


(41) 


where 


°ak  '  VI  *2  \  ■ 


(42) 
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In  (41),  we  have  0  <  a  <  K  and  it  is  understood  that  ct^q  =  0  whenever 
|m|  >  N.  Ogj-  is  readily  found  from  (23)  and  (28): 


(43) 


(44) 


B  is  a  vector  obtained  from  the  right  hand  side  of  (41),  and  M  is  the 
matrix  of  coefficients  obtained  from  the  left  side  of  (41).  We  give  an 
example  of  M  and  B  below  in  Section  III.  The  numerical  results 
presented  below  have  been  obtained  by  solving  (41)  numerically. 

For  this  method  to  yield  accurate  results,  N  must  be  chosen  large 
enough  that  at  least  2LJ  of  the  are  essentially  zero  on  each  end  of 

the  vector  A  in  (46).  If  this  vector  has  2LJ  zeros  on  each  end,  the 
solution  obtained  from  (41)  will  yield  the  same  result  as  the  solution 
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of  the  infinite  system  in  (38).  In  practice,  a  suitable  value  of  N  may 
be  determined  by  increasing  N  until  the  A  vector  has  2LJ  negligible 
terms  on  each  end  and  until  the  values  of  the  ap^  in  the  middle  of  the 
vector  A  are  unaffected  by  further  increases  in  N.  Experience  shows  how 
large  N  must  be  in  specific  cases*.  Once  the  ap  ^  have  been  found,  the 
Cn  may  be  found  from  (40)  and  W(t')  from  (39). 

The  time-varying  array  weights  have  two  effects  on  array 
performance.  They  cause  the  array  to  modulate  the  desired  signal,  and 
they  cause  the  array  output  SINR  to  vary  with  time.  In  the  remainder  of 
this  section,  we  discuss  the  calculation  of  these  effects. 

Given  a  time-varying  weight  vector  W(t'),  the  desired  signal 
component  of  the  array  output  is 


s~d(f)  =  Ad  wT(t')  Ud 


+  V 


(47) 


1  2  ~ 
where  <dq  =  oiQ/ko  .  To  study  the  modulation  on  sd(t),  we  define 

ad(f)  eJnd(t,)  =  A/(f)Ud  .  (48) 

Then  ad(t‘)  =  Ad|WT(t')Ud|  is  the  envelope  modulation  and  nd(t‘)  s 
«wT(t')Ud  is  the  phase  modulation.  Furthermore,  we  define  adn(t’)  to  be 


♦The  numerical  solution  of  (41)  was  done  using  a  Gauss  elimination 
routine  with  full  pivoting  [9],  We  have  not  examined  the  eigenvalue 
behavior  of  M  in  detail,  but  we  did  not  experience  any  convergence 
problems  or  other  unusual  behavior  due  to  ill -conditioning  with  this 
approach. 
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the  envelope  normalized  to  its  value  in  the  absence  of  interference, 
i  .e.. 


= 


ad(t‘) 
Ad ' WoTud i 


(49) 


where  W0  is  the  steady-state  weight  vector  that  would  occur  without 
i nterference, 

W0  =  (^d  +  <&n)-1S  •  (5°) 


(4(1,  $n  and  S  are  given  in  (17),  (20)  and  (22).)  Our  results  below  are 
presented  in  terms  of  aCjn(t'),  because  the  effect  of  the  interference 
can  be  seen  directly  by  comparing  adn(t')  with  unity. 

The  output  desired  signal  power  is 

Pd(t,)  =7E[lsd(t'),2l  =?Ad  lwT(t,)  Ud*2  *  (51) 

The  output  interference  power  is 

Pi  (t1)  =  |  A2  |WT(t')  ui (t  * ) ! 2  ’  (52) 

where  U-j  ( t ‘ )  is  the  vector  of  phasors  given  in  (11)  but  written 

in  terms  of  normalized  time  t1  =  ka^t.  The  output  thermal  noise  power 

is 


Pn(f)  =  WT(t’)  W*(t') 


(53) 
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From  these  quantities,  the  output  interference-to-noise  ratio  (INR), 
INR  =  pi (*• '  )  .  (54) 

WT 

and  the  output  signal -to-interference-plus-noise  ratio  (SINR), 


may  be  computed  as  functions  of  t'. 

In  the  next  section,  we  present  an  example  using  this  technique  and 
discuss  the  effects  of  a  phase  modulated  interference  signal  on  the 
adaptive  array. 


III.  AN  EXAMPLE 

Consider  an  array  with  two  elements,  so  J=Z  and  K=l.  Let  the 
interference  have  sinusoidal  phase  modulation, 

Y.  (t1)  =  Bsinoit  ♦  (56) 

1  m 


where  6  is  the  maximum  phase  deviation  (or  the  modulation  index)  and  u' 

m 

is  the  normalized  modulation  frequency.  Then  one  finds  from  (19)  and 
standard  trigonometric  identities  that 
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where 


8  =  28sin 


oi't! 
m  i 


(and  T.!  =  kon).  Using  the  Fourier  Series  expansion  [101, 


jzsinp  jtp 

e  =  l  0j(z)e 

£=-<*> 


we  find 


?  e 
i 


jx^ft) 


r  j?V  t1 

If  e  m 

J,  =  _oo  t01 


with 

f,0i  ■  wn  e'jl  • 

From  (34),  we  also  have 

fuo  =  f-*oi  =  f*01  ’ 


since 

=  M)*  Jt(B')  . 

To  truncate  the  series  in  (60),  we  use  the  fact  that 
Jj(z)  5  0  for  | t|  >  z  +  1 


(58) 


(50) 


(60) 


(61) 


(62) 


(63) 


(64) 
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Thus,  we  approximate 


jxoi(t')  L  jtu't' 

5.e  =  l  f  e  m 

i  «,=-!_  W1 


(a)1  T* 

L  =  {g*  +  1}  =  {28  sin  _!!L-1  +  1} 


where  (r)  denotes  the  smallest  integer  greater  than  r. 

Writing  out  (41)  with  L=2  as  an  example  and  with  K=1  gives  the 
following  matrix  of  coefficients  M  in  (45): 


..000  F_2  0  F_!  d2  Z  0  Fj  0  F2 

;....0F*  0  0  Z*d2  F^o  f*2  0  0  0 . 

....  000  F_2  0  F_i  dj  Z  OF!  0  F2  0  .... 

....  0  F*  0  F*  0  Z*  di  f*  0F*  0  0  0  .... 

2  1  -1  -c 

....  000  F_2  0  F_i  do  Z  0  Fj  0  F2  0  .... 

....  0  F*  0  F*  0  Z*  d0  F*j  0  f*_20  0  0  .... 

_  00  0  F_2  0  F_id_i  Z  0  Fi  0  F2  0  .... 

_  0  F*  0  F*  0  z*  d-l  F*j  0  F*2  0  0  0 
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The  vectors  A  and  S  in  (45)  are 


“N.O 


a2,0 

C',1  j 

al,0 

al,l 

“O.O 

<*0,1 

a-l,0 

a-l,l 


where 

f£  =  tOl  » 

which,  from  (34),  also  implies 
★  * 

F-4  =  f-Sl01  =  ftlO  ’ 
and  where 


dn  =  f000 


+  Q00  +  jn“m  =  1  +  5d 
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and 


Z  =  Q, 


01 


001 


=  C, 


j  ( 4>i  -  *d ) 


Vo(S') 


(71) 


To  obtain  the  results  presented  below,  we  have  solved  this  system  (with 
the  appropriate  values  of  L)  numerically  for  the  a^,  and  then 
evaluated  Cn  and  W(t')  from  (40)  and  (39). 

Let  us  now  discuss  the  results  obtained  from  these  calculations . 

In  part  A  below,  we  show  typical  curves  of  desired  signal  modulation, 
output  INR  and  SINR  as  functions  of  time.  In  parts  B-E,  we  describe  the 
effect  of  each  signal  parameter  on  the  desired  signal  modulation.  In 
part  F,  we  assume  the  array  is  used  in  a  DPSK  (differential  phase  shift 
keyed)  digital  communi cati on  system  [111  and  show  how  the  bit  error 
probability  is  affected  by  the  interference  signal  parameters  C-j ,  B  and 


A.  Typical  Waveforms 

Figure  2  shows  typical  curves  of  the  normalized  envelope  modulation 

u)'  9 

adn (t  * )  for  9d=30°,  9.=45°,  5d=10  dB,  ^=30  dB,  f^  =  ^  =  10  , 

f '  =  JL  -  10®  ,  and  6=10^ ,  10®,  104  and  10®.  These  curves  illustrate 
o  2? 

that  the  envelope  modulation  on  the  desired  signal  can  be  substantial 
for  larger  values  of  6. 

It  turns  out  that  for  this  2-element  array,  angle  modulated 
interference  does  not  cause  phase  modulation  on  the  desired  signal. 

This  result  was  discovered  by  calculating  n<j ( t ' ) ,  for  numerous  values 
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j3  =  I05 

— | - 1 - 1 - 1 - r— i - 1 - 1 - 1 - 1 - 1 - 1 - 1 - r 

0.0025  0.005  0.0075 

TIME 

Figure  2.  adn{t')  versus  time. 

9d=30°,  01=45°,  f-d  =  10  dB,  =30  dB 


of  9.,  6,  and  ? . .  These  calculations  show  that  ( t 1 )  does  not 
change  with  time  regardless  of  the  signal  parameters.  The  reason  for 
this  behavior  is  as  follows.  Roth  the  amplitude  and  phase  of  each 
weight  in  the  array  vary  sinusoidally  with  time.  The  amplitudes  of  the 
two  weights  are  equal  at  each  instant  of  time.  The  phase  angle  of  each 
weight  contains  a  term  constant  with  time,  which  depends  on  interference 
arrival  angle,  and  a  term  that  varies  sinusoidally  with  time.  The 
sinusoidally  varying  term  is  180°  out  of  phase  on  the  two  weights.  When 
a  desired  signal  is  passed  through  these  weights  the  phase  modulation 
produced  on  the  desired  signal  by  one  weight  is  180°  out  of  phase  with 
that  produced  by  the  other  weight,  and  the  sum  signal  at  the  array 
output  contains  no  phase  modulation. 

Figures  3  and  4  show  typical  curves  of  the  output  INR  and  SINR  as 
functions  of  time,  over  one  period  and  for  the  same  signal  parameters  as 
in  Figure  2.  Figure  3  shows  that  the  average  INR  increases  and  the 
array  SINR  decreases  as  8  increases.  Also,  as  8  increases,  the  INR  and 
SINR  variations  with  time  are  more  pronounced.  At  low  8,  the  INR  and 
SINR  are  constant,  because  the  array  feedback  tracks  the  incoming  phase 
modulation. 

In  general  one  finds  that  the  array  behavior  changes  substantially 

as  the  signal  parameters  8  0^ ,  ,  f^  and  8  are  varied.  In  parts 

8-E  of  this  section,  we  examine  the  effect  of  each  signal  parameter  on 
the  desired  signal  modulation.  In  part  F,  we  show  how  bit  error 
probability  is  affected  by  the  angle  modulated  interference  when  the 
array  is  used  in  a  DPSK  communication  system. 


To  describe  the  desired  signal  modulation,  we  shall  define  three 
quantities.  First,  we  let  amin  and  amax  be  the  minimum  and  maximum 
values  of  aq^t1)  during  the  modulation  period.  Then,  we  define 

m  =  amax  ~  amin  .  (72) 

amax 

m  is  the  envelope  variation  normalized  to  its  peak.  It  may  be  thought 
of  as  "fractional  modulation",  analogous  to  percentage  modulation  in  AM. 
In  Parts  B-E  below,  we  describe  the  effect  of  each  signal  parameter  on 
amax  an^ 

B.  The  Effect  of  Angle  of  Arrival 

Desired  signal  modulation  effects  are  small  when  ^  is  far  from  9d. 
When  9-j  is  close  to  0q,  the  envelope  variation  m  is  large  and  the  peak 
amax  is  small.  However,  when  9-j  is  equal  to  9q,  the  desired  signal  is 
nulled,  so  m  drops  to  a  small  value. 

Figures  5  and  6  show  typical  curves  of  m  and  amax  as  functions  of 

9.  for  9d=0°,  Cd=10  dB,  C.=40  dB,  f^=103  and  f^lO8.  Four  different 
curves  are  shown,  for  8  =  2x10^,  4x10^,  8x10^  and  10^.  It  is  seen  that 
m  is  large  and  amax  is  small  when  9}  is  near  9d.  We  note  also  that  when 
9-j  is  extremely  close  to  6d,  m  drops  to  zero.  This  behavior  occurs 
because  as  9j  becomes  very  close  to  9d,  the  peak-to-peak  variation 
(amax~amin)  approaches  zero  more  quickly  than  amax  does,  so  m  goes  to 
zero.  However,  this  behavior  of  m  is  no  importance,  since  the  desired 
signal  is  nulled  anyway. 
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/  z: 


9d=0°.  5d=10  dB,  Ci=40  dB 


9.  ( DEGREES ) 

i 

Figure  6.  amax  versus  9^. 

9<i=0°.  ^=10  dB,  e,  =40  dB 


C.  The  Effect  of  Modulation  Index  and  Frequency 

The  peak  a  _  is  large  at  low  f'  and  drops  as  f*  increases.  The 
r  ma  x  3  m  K  m 

higher  B,  the  farther  a  drops  at  high  f'.  The  variation  m  is  small 

max  m 

at  low  f',  peaks  at  intermediate  values  of  f'  and  drops  as  f'  increases, 
m  r  m  m 

For  low  and  intermediate  f ’ »  m  is  largest  for  high  B.  At  high  f',  m 

m  m 

drops  more  rapidly  with  larger  values  of  B. 

Figures  7  and  8  illustrate  these  results  for  the  case  9^=30°, 

9^=45°,  ?d=10  dB,  £^=30  dB,  ^=10®  and  Tor  values  of  B  between  2x10^ 

4  1 

and  i0  .  It  is  seen  in  Figure  7  that  a  is  large  for  low  f  and 

drops  as  f^  increases.  This  behavior  is  due  to  the  array  speed  of 

response,  which  limits  its  ability  to  track  the  modulation  at  high  f\ 

m,  seen  in  Figure  8,  peaks  at  intermediate  f\ 

D.  The  Effect  of  Interference-To-Noi se  Ratio 

For  all  values  of  f^,  the  variation  m  is  largest  for  high  I NR . 

a  is  hiqh  at  low  f'  and  drops  at  high  f'. 
max3  m  r  3  m 

I 

Figures  9  and  10  show  the  behavior  of  a  ,  and  m  versus  f  for 

max  m 

9d=30°,  9. =45°,  Cd=10  dB,  S=104,  f^=108  and  for  four  values  of  ^ 

between  20  and  45  dB.  It  may  be  seen  in  Figure  9  that  a  peaks 

max 

slightly  at  high  INR,  for  values  of  f^  around  400.  This  behavior  is  due 
to  the  change  in  array  speed  of  response  with  INR. 
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m 


a 

c\j 


Figure  8.  m  versus  f ' . 

m 

9d=30°,  6-j  =45°,  £d=10  dB,  £1=30  dB, 

f;-io8  . 
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E.  The  Effect  of  Desired  Signal -To-Noise  Ratio 


m  is  largest  and  amax  is  smallest  for  low  5^.  As  is  increased, 
m  decreases  and  amax  increases. 

Figures  11  and  12  illustrate  this  behavior.  Figure  11  shows  m  and 

Figure  12  shows  a  ,  both  versus  f'  for  9  =30°,  8. =45°,  £.=30  dB, 
max  m  d  i  i 

f p=10^,  8=lo\  and  for  four  values  of  between  10  and  40  dB.  It  is 

seen  that  m  peaks  at  intermediate  f'  with  the  value  of  f*  at  the  peak 

m  m 

dependent  on  . 

F.  Bit  Error  Probability 

To  evaluate  the  effect  of  the  time-varying  SINR,  we  have  computed 
the  bit  error  probability  when  the  desired  signal  is  DPSK  biphase 
modulated  signal.  We  assume  the  desired  signal  bit  duration  is  large 
compared  with  the  propagation  time  across  the  array  but  small  compared 
to  the  interference  modulation  period  Also,  we  assume  the 

reference  signal  is  replica  of  the  desired  signal.  Under  these 
assumptions  the  array  weights  with  the  DPSK  biphase  modulated  desired 
signal  are  the  same  as  those  produced  with  a  CW  signal.  The  effective 
bit  error  probability  Ve  has  by  obtained  by  averaging  the  instantaneous 
bit  error  probability  over  one  period  o*  the  interference  modulation. 
The  instantaneous  bit  error  probability  Pe ( t ' )  is  related  to  the  SINR 
by  [12], 
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(73) 


?e(f) 


=  1  e-SINR (t ' ) 


Figure  13  shows  typical  curves  of  Pg  as  a  function  of  f^  for 

9^=30°,  9. =45°,  5^=10  dB,  £.=30  dB,  f0=lO8  and  for  four  values  of  8 
3  4  1 

between  2x10  and  10  .  At  low  f  ,  where  the  interference  bandwidth 

m 

is  too  small  to  cause  any  degradation  in  the  array  performance,  stays 
constant  at  a  value  that  does  not  depend  on  6.  This  value  of  is 
determined  simply  by  the  input  I NR .  As  f^  increases  the  interference 
bandwidth  has  more  effect.  At  high  f^,  Pg  approaches  0.5. 

Figure  14  shows  typical  curves  of  P0  versus  f^  for  different  values 
of  £.  and  for  9^=30°,  9^=45°,  £d=10  dB,  f0=lO8  and  8=104.  These  curves 
illustrate  the  dependence  of  P0  on  the  interference  power.  At  low  f', 
the  interference  has  small  bandwidth  and  depends  slightly  on  Si . 

This  behavior  is  the  same  as  with  CW  interference. 


IV.  CONCLUSIONS 

We  have  developed  a  method  for  finding  the  periodic  steady-state 
weights  of  an  adaptive  array  receiving  an  interference  signal  with 
periodic  angle  modulation.  We  then  used  this  method  to  study  the 
behavior  of  a  two-element  array  subjected  to  an  interference  signal  with 
sinusoidal  phase  modulation.  This  interference  causes  the  array  weights 
to  cycle  sinusoidally  with  time,  both  in  magnitude  and  phase.  The 
modulated  weights  cause  the  desired  signal  to  have  envelope  (but  not 
phase)  modulation  at  the  array  output. 
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Figure  13.  Bit  error  probability  versus  f'. 

m 

9d=30\  61=45°,  ed-10  dB, 


£. =30  dB,  f^=10S  • 
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The  envelope  modulation  is  largest  when  the  interference  arrives 
close  to  the  desired  signal  in  space.  The  envelope  modulation  increases 
with  the  bandwidth  of  the  input  interference  signal.  Phase  modulated 


interference  also 
the  array  is  used 
increases  the  bit 


causes  the  array  output  SINR  to  vary  with  time.  When 
in  a  digital  communication  system,  such  SINR  variation 
error  probability. 
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